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Integrable systems with boundary interactions are one of the recent achievements which the authors think deserve 
careful investigations. Here we only concern with 1-dimensional open-boundary lattice integrable models of strongly 
correlated electrons. Such lattice models can be treated by Sklyanin's boundary quantum inverse scattering method 
(QISM) or its generalizations More specifically, we present integrable boundary terms for the eight-state 

version of the supersymmetric U model recently introduced in Q . The bulk Hamiltonian describes a supersymmetric 
electron model with correlated single-particle and pair hoppings as well as uncorrelated triple-particle hopping. So 
the model on an open lattice, which we consider here, involves many physically interesting processes with boundary 
interactions. The boundary model is solved by means of the coordinate Bethe ansatz method and the Bethe ansatz 
equations are derived. 

Let Cj a (cj, a ) denote fermionic creation (annihilation) operator which creates (annihilates) an electron of species 

a (= +,0,—) at site j. They satisfy the anti-commutation relations given by {cJ„,Cj,^} = 5ijS a p, where i,j — 
1, 2, ■ • • , L and a,/3 = +, 0, — . We consider the following Hamiltonian with boundary terms 

L-l 

H = j2 #j™v k i + i^ oundary + H h R oundaiy , (i) 

i=i 

where -ffj^+i is the Hamiltonian density of the eight-state supersymmetric U model || 

H jf+i(9) = ~ X^L ^ 1 -" + h ' c ') ex P { ~\ + n i+ 1 'P) + \ H ( n j,P n 3,i + n j+ i,f3n j+1 ^) 



2{i+V) ^ ( C L4/3 C J + 1 .' 3C J + 1 ^ +h - C -) eX p{-f( n J,7+ n j + l,7)j 

• (g + l)(g + 2) (4 + c ioi-^ +1 ,_c J+1 , oCj+1 , + +h.c.) 
^2(nj, a +n j+ha ) - (n jta n jt p + n j+ha n j+hl3 ) 

2 



(g+l)(g + 2) 



• "j.ullj. + n j+l, + n j+l,Q n ] + l-) (2) 



with 



In -Jr. C-ln(.9 + l)-iln.g( 5 + 2), £ = -ln ^ (3) 
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j rr boundarv rr boundary i 1 , 

and H L , H R are boundary terms 

rjboundary ^9 ( ^ / , , \ . ^ 

/f £ • = - - _ ^ I ni - — (m+riio + ni ni_ + nx+rix-) + ^ - - ^ n 1+ n 1Q n^ 

^boundary _ — i?_ _ JL( nL+n , i0 + nLO n L _ + n L +n L -) + . 8 r n L+ n L0 n L - ) . (4) 
* _ 4+ \ 4+ 4+t/ + 4+J / 

In the above equations, nj Q is the number density operator n JQ = cj Q c JQ , nj = nj_|_ + Ujo + Tij—j £± are some 
parameters describing the boundary effects. 

As was shown in [||, the supersymmetry algebra underlying the bulk model is gl(3\l). The boundary terms may 
spoil this <?^(3|1) supersymmetry, leaving the boundary model with a smaller symmetry algebra. If one projects out 
any one species, then the projected Hamiltonian is nothing but that of the supersymmetric U model with boundary 
terms Q with the following identification of the parameters in the two models: U = i^qrp 

We will establish the quantum integrability of the Hamiltonian for the boundary eight-state supersymmetric U 
model ([!]) by showing that it can be derived from the (graded) boundary quantum inverse scattering method. Let 
us first of all recall that the Hamiltonian of the eight-state supersymmetric U model with the periodic boundary 
conditions commutes with the bulk transfer matrix, which is the supertrace of the monodromy matrix T(u), 

T(u) = R 0L (u)---R 01 (u). (5) 

where the subscript denotes the auxiliary superspace V — C 4 ' 4 . It should be noted that the supertrace is carried 
out for the auxiliary superspace V. The R-matrix R(u) = PR(u), where P is graded permutation operator, is given 
by 

( > u-2g (u-2g){u-2g-2) {u - 2g)(u - 2g - 2)(u - 2g - 4) { > 

where Pk, k — 1,2,3,4, are four projection operators whose explicit formulae can be found in Q. The elements of 
the supermatrix T(u) are the generators of an associative superalgebra A defined by the relations 

Rx 2 (u x - u 2 ) T (mi) T (u 2 ) =T (u 2 ) T (ux)Rx2(ux ~ u 2 ), (7) 
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where X= X ® 1, X= 1 <&X for any supermatrix X S End(V). For later use, we list some useful properties enjoyed 
by the R-matrix: (i) Unitarity: Rx 2 {u)R 2 x{— u) = 1 an d (ii) Crossing-unitarity: R s ^ 2 2 (—u + 4)i?2* 1 (u) — p{u) with 
p(u) being a scalar function. 

In order to construct intcgrablc electronic models with open boundary conditions, we introduce the following graded 
reflection equations (REs) that the so-called boundary K-(super)matrices K±(u) satisfy 

1 2 2 1 

#12 (wi - u 2 ) K- {ui)R 21 (ux + u 2 ) K- (u 2 ) =K- (u 2 )Ri 2 (ui + u 2 ) K- (ui)R 21 (ux - u 2 ), (8) 



R s 21 llst2 (-ux+u 2 ) 4* 1 (u^Rui-u.-^+^K^ (u 2 ) 

=K\ t2 (u 2 )R 21 (-ux-u 2 + 4)K^ (ux)R s xl list2 (-ux+u 2 ), (9) 

where the supertransposition st^ (/x = 1, 2) is only carried out in the /i-th factor superspace of V <S) V , whereas ist^ 
denotes the inverse operation of st^. Following Sklyanin's arguments |EJ, one may show that the quantity 71 (it) given 
by 

71(u) = T(u)K-(u)T-\-u) (10) 

satisfies the same relation as K—(u): 

1 2 2 1 

R 12 (ux - u 2 ) T- (u x )R 2 x{ux + u 2 ) T- (u 2 ) =T_ {u 2 )R 12 (ux + u 2 ) T- (ux)R 2 x{ux ~ u 2 ). (11) 
Thus if one defines the boundary transfer matrix t(u) as 
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t(u) = str{K+(u)T-{u)) = str (K + (u)T{u)K_{u)T- l (-u)) , 
then it can be shown j|] that 

[r(ui),r(«2)] = 0. 



(12) 



(13) 



We now solve (||) and d) for K+(u) and K_{ u). Let us restrict ourselves to the diagonal solutions. Then, one may 
check that the matrix K_(u) given by 



£_(2-£_)(2 + £_) 



/ A_(u) 
























B_(u) 












B_(«) 












C_(«) 












CL(«) 












C-(u) 




(14) 



£>_(«) / 



where 







= (■ 




+ «)(2 




+ u)(-2 




+«), 




(«) 


= (- 




+ «)(2 




-u)(-2 




+«), 


C_ 


(«) 


= (■ 




-u)(2 




-u)(-2 




+«), 




■(«) 


= (- 




-«)(2 




-«)(-2 




- u), 



(15) 

satisfies (|8|). The matrix K+{u) can be obtained from the isomorphism of the two REs. Indeed, given a solution 
K-(u) of (pi), then K + (u) defined by 



K s Hu) = K_(-u + 2) 



is a solution of (O). The proof follows from some algebraic computations upon substituting (16) into 



(16) 

and making 



use of the properties of the R-matrix . Therefore, one may choose the boundary matrix K + (u) as 



K+{u) 



( A+(u) 












B+(u) 












B+(u) 












B+(u) 












C+{u) 












C+(u) 












C+(u) 












D+{u) j 



(17) 



where 



A+(u) = (-2 fl + 2 + £+ - u)(-2g + £+ - u){-2g - 2 + £+ - u), 
B+{u) = {-2g - 2 + + u){-2g + £+ - u)(-2g - 2 + - u), 
C+{u) = {-2g - 2 + £+ + u){-2g - 4 + £+ + u){-2g -2 + £+-u), 
D+{u) = (-2g-2 + Z + +u)(-2g-4: + Z + +u)(-2g-6 + Z + +u). 



(18) 



Now it can be shown that Hamiltonian (|l|) is related to the second derivative of the boundary transfer matrix t(u) 
(up to an unimportant additive constant) 



H = 2g H u , 
r"(0) 



L-l 



H R = 



A{V + 2W) 



2 



(0) 



1 



2(V + 2W) 



str K+ (0)G LO 



+2 str K'+ (0)Ff + str if + (0) (ff*, 



(19) 
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where 



V = str K' + (Q), W = str (Jc+ (0)H* 
H% = PijRij (0), Gij = PijU'lj (0). (20) 

Here Pjj denotes the graded permutation operator acting on the i-th and j-th quantum spaces. (|l9|) implies that the 
boundary eight-state supersymmetric U model admits an infinite number of conserved currents which are in involution 
with each other, thus assuring its integrability. ft should be emphasized that Hamiltonian ([[]) appears as the second 
derivative of the boundary transfer matrix t(u) with respect to the spectral parameter u at u — 0. This is due to 
the fact that the supertrace of ^+(0) equals to zero. The reason for the zero supertrace of K + (0) is related to the 
fact that the quantum space is the 8-dimensional typical irreducible representation of gl(3\l). A similar situation has 
already appeared in many spin chain models ■ 

Having established the quantum integrability of the boundary model, we now solve it by using the coordinate space 
Bethe ansatz method. The whole procedure is similar to that for other models ||f7| . The Bethe ansatz equations are 



e 



a=i 0j - Air' - *c/2 9j + \£ } - ic/2 



" X W _ g . + lc/2 A (D + g. + ic/2 Ml A (D _ A W + lc Ag + A W + lc ^ A W - Ag - *c/2 A< 1} + \f ] - ic/2 
/i A« - 6, - *c/2 A« + 6, lc /2 A« - A< 1} - zc A« + A« - zc ^ A« - A< 2) + zc/2 A« + A< 2) + zc/2 ' 

a = !,-■■, Mi, 

^ Ag _ A W + ic/2 Ag + Ag + zc/2 _ g Ag - A< 2) + zc A< 2) + A< 2) + zc 

7 = !,-••, M a , (21) 



where 



CO;p) = - — > pi = -i~o — px, = -1-75 — 7-, (22) 

pe lK 2 — £_ 2 — ^ 



and c = e v — 1; the charge rapidities 9j = 6(kj) are related to the single-particle quasi-momenta kj by 9(k) = ~ tan(-|). 
The energy eigenvalue E of the model is given by E = —2 cos kj (modular an unimportant additive constant 

coming from the chemical potential term). 

In conclusion, we have studied integrable open-boundary conditions for the eight-state supersymmetric U model. 
The quantum integrability of the system follows from the fact that the Hamiltonian may be embedded into a one- 
parameter family of commuting transfer matrices. Moreover, the Bethe ansatz equations are derived by use of the 
coordinate space Bethe ansatz approach. This provides us with a basis for computing the finite size corrections (see, 
e.g. HQ) to the low- lying energies in the system, which in turn allow us to use the boundary conformal field theory 
technique to study the critical properties of the boundary. The details will be treated in a separate publication. 
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